Draft version February 1, 2008 

Preprint typeset using IiT^C style emulateapj v. 04/03/99 



COSMOLOGICAL EFFECTS OF THE FIRST STARS: EVOLVING SPECTRA OF POPULATION III 

JASON TUMLINSON 1 , J. MICHAEL SHULL 2 , & APARNA VENKATESAN 3 
Center for Astrophysics and Space Astronomy, Department of Astrophysical and Planetary Sciences, 
University of Colorado, Boulder, CO 80309 
Draft version February 1, 2008 

ABSTRACT 

The first stars hold intrinsic interest for their uniqueness and for their potentially important contributions to 
galaxy formation, chemical enrichment, and feedback on the intergalactic medium (IGM). Although the sources 
of cosmological reionization are unknown at present, the declining population of large bright quasars at redshifts 
z > 3 implies that stars are the leading candidates for the sources that reionized the hydrogen in the IGM by 
z ~ 6. The metal-free composition of the first stars restricts the stellar energy source to proton-proton burning 
rather than the more efficient CNO cycle. Consequently they are hotter, smaller, and have harder spectra than 
their present-day counterparts of finite metallicity. We present new results from a continuing study of metal-free 
stars from a cosmological point of view. We have calculated evolving spectra of Pop III clusters, derived from 
a grid of zero-metallicity stellar evolutionary tracks. We find that H-ionizing photon production from metal-free 
stellar clusters takes twice as long as that of Pop II to decline to 1/10 its peak value. In addition, metal-free stars 
produce substantially more photons than Pop II in the He II (E > 4 Ryd) continuum. We suggest that large Lya 
equivalent widths (WLy Q > 400 A) may provide a means of detecting metal-free stellar populations at high redshift, 
and that He II recombination lines (A1640, A4686) may confirm identifications of Population III. While Pop III 
clusters are intrinsically bluer than their Pop II counterparts, nebular continuum emission makes up this difference 
and may confuse attempts to discern Pop III stars with broadband colors. In a companion paper, we explore the 
consequences of evolving spectra of Pop III for the reionization of the IGM in both H and He. 



1. INTRODUCTION 

The first stars and their role in galaxy formation and the 
evolution of the intergalactic medium (IGM) define one of the 
frontiers of modern cosmology. Yet we do not know when 
the first stars formed, where they now reside, what happened 
to them over time, or even where to look for them. The Big 
Bang paradigm requires that the first stars were composed of 
H and He with only trace light elements and negligible carbon 
(Schramm & Turner 1998). Metal-free stars have unique prop- 
erties owing to this primordial composition: they are hotter and 
smaller than their metal-enriched counterparts (Tumlinson & 
Shull 2000, hereafter TS). These features have interesting im- 
plications for their formation, their environment, and particu- 
larly for their spectra. Some key open issues are the formation 
epoch of the first stars (Ricotti, Gnedin, & Shull 2002), the mass 
distribution for stars formed from primeval gas (Bromm, Coppi, 
& Larson 2001), and their contribution to the reionization of H 
and He in the IGM. The recent suggestions by Lanzetta et al. 
(2002) that star formation may become more widespread in the 
distant past, the discovery of a gravitationally-lensed galaxy at 
Z = 6.56 (Hu et al. 2002), and the apparent completion of H 
reionization by z — 6 (Becker et al. 2001; Djorgovski et al. 
2001) and He at z — 3 (Kriss et al. 2001) make detailed exami- 
nation of the first stars particularly timely. 

With our present tools we can see back to z ~ 6, roughly 1 
Gyr after the Big Bang (a universe with il m = 0.3, ft\ = 0.7, 
and Ho = 70 km s" 1 Mpc" 1 is assumed here). Imaging and 
spectroscopic surveys that sample wide swaths of sky for tell- 
tale emission lines (Lya emitters; Hu et al. 1999; Rhoads et 
al. 2000) or peculiar broadband colors (Lyman-break galaxies, 



Steidel et al. 1999) have discovered more than 1000 galaxies at 
z ~ 3. Searches such as the Sloan Digital Sky Survey (SDSS) 
are expanding the list of known QSOs at z > 5 and enabling 
high-redshift probes of the IGM (Fan et al. 2001; Schneider 
et al. 2002). Ground-based spectra can determine the star for- 
mation rate, metallicity, and mass of galaxies at z = 3 (Pettini 
et al. 2002). The first stars are believed to have formed be- 
tween z = 20-30 (Gnedin & Ostriker 1997; Ricotti et al. 2002), 
while galaxies have been discovered up to z ~ 6. The next ob- 
servational frontier is the direct observation of the first stellar 
generation. 

This work is part of an ongoing effort to examine the first 
stars from a unified stellar and cosmological point of view. We 
began with zero-age models (TS) of metal-free stars that re- 
vealed their potential importance to H and He reionization in 
the IGM. We followed that with a detailed calculation of the 
He II recombination line fluxes from Population III clusters, 
assessing the prospects for discovering this emission as a dis- 
tinctive signature of the first stars (Tumlinson, Giroux, & Shull 
2001, hereafter TGS). We are concerned here with the proper- 
ties of the first stars that affect their observational signatures 
and feedback on the IGM. These effects depend on the answer 
to a fundamental question about the first stars: what is the time 
evolution of their structure and radiation? We update the work 
of TS to compute new evolving spectra of metal-free stars and 
to assess the effects of evolution on the emergent spectra of 
Population III clusters. 

In this paper, we use the terms "Population III" and "metal- 
free" interchangeably. We also note that there are two distinct 
usages of the term "Pop III" that have appeared in the litera- 
ture. One usage connotes small protogalactic clusters that have 
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Mum — 10 Mq and that must cool with molecular hydrogen 
lines in the absence of metals. We adhere to the stellar def- 
inition, in which "Pop III" and metal-free stars are identical 
regardless of location. 

In § 2 we discuss the major components of our stellar evolu- 
tion calculations. In § 3 we present the stellar evolution tracks 
that enable our evolving spectra and observational predictions. 
In § 4 we describe the modeling of the stellar atmospheres and 
construct evolving spectra for synthetic Pop III clusters. In § 5 
we make detailed predictions of the observational signatures 
of metal-free stars and describe the realistic prospects for their 
discovery. In § 6 we draw some general conclusions. Paper 
II (Venkatesan, Tumlinson, & Shull 2002) presents the results 
of a cosmological reionization model that incorporates Pop III 
spectra as ionizing sources. 

2. STELLAR EVOLUTION AT ZERO METALLICITY 

2.1. Method 

To study the evolution of Population III stars, we must con- 
struct realistic numerical models that account for the physical 
processes that are thought to be relevant in real stars. The ba- 
sic theory of stellar structure is quite well developed and has 
been tested against observations for decades. Occasional con- 
flicts with detailed observations lead to refinements of the the- 
ory, usually in the form of additional input physics (convective 
overshoot, mass loss, and rotation, for example). Of course, 
no such fine calibration of models for Pop III stars is yet pos- 
sible, since we have not discovered any of these stars or any- 
thing believed to be their remnants. We must make the basic as- 
sumption that the framework of stellar evolution theory applies 
equally well to the unusual regime of zero metallicity. Under 
this assumption, our Population III models incorporate the same 
paradigmatic assumptions and physical ingredients that under- 
lie all modern stellar evolution codes. We must note here, how- 
ever, that since our ultimate goal is a calculation of the evolving 
spectra, we are not concerned with internal processes, particu- 
larly late in the stellar lifetime, that do not affect the spectral 
evolution. Those readers who are interested in the details of 
the later stages of Pop III evolution are referred to the tracks 
by Marigo et al. (2000) and Siess, Livio, & Lattanzio (2002). 
Schaerer (2002) has compiled tracks from a variety of sources 
(El Eid et al. 1983; Marigo et al. 2000; Klapp 1983) in an in- 
vestigation into the properties of metal-free stars that parallels 
some of our own results. Finally, we note that the following ex- 
position of our method is included for readers interested in the 
details. The basic results of our study can be found by skipping 
to § 3. 

In our code we solve the fundamental stellar structure equa- 
tions: 
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where radius R, luminosity L r , pressure P, and temperature T 
are the "primary variables" and the subscript r implies that the 
variable represents a cumulative quantity up to radius R. These 



equations are transformed to a Lagrangian mass grid that does 
not change with total stellar mass. The mass density p, opac- 
ity k, energy generation rate e, thermodynamic derivatives V ra d 
and V a d, and the specific entropy S are "secondary variables" 
which are themselves functions of T, P, L, R, and the chemi- 
cal abundances by number, denoted by X t . For radiative energy 
transport, we use: V = V rac j = (dlnT /dlnP)^ and in convec- 
tive regions we use V = (<f lnT/a! lnP) ad . The adiabatic deriva- 
tive is calculated for each mass point during the evaluation of 
the equation of state. 

We solve these coupled equations with a variant of the relax- 
ation method described definitively by Kippenhahn, Weigert, & 
Hofmeister (1967) and still in extensive use today. This method 
transforms the main differential equations into difference equa- 
tions by expanding them in a Taylor series about small correc- 
tions to the initial guesses for the primary variables. The cor- 
rections are derived iteratively by a routine adapted from Press 
etal. (1987; Chapter 16.3). 

2.2. Physical Ingredients of the Models 

2.2.1. Equation of State 

The equation of state (EOS) is a numerical function that re- 
lates the primary variables pressure P, temperature T, and the 
secondary variable density p. The EOS is that of an ideal 
monatomic gas with radiation pressure. Partial ionization of 
the metallic species is included but is generally unimportant at 
T > 10 7 K in the stellar core. 

2.2.2. Opacity 

Our code computes the Rosseland mean opacity k = 
n(p,T,X,Z) from the OPAL tables produced at the Lawrence 
Livermore National Laboratory (Iglesias & Rogers 1992). The 
code calls a subroutine provided publicly by the OPAL project 
to interpolate in the tables and derive k for a specific combina- 
tion of parameters (p, T,X,Z). The tables are provided at fixed 
grid points in T and R (R = p/(T/\0 6 K) 3 , not to be confused 
with the radius) and there are multiple tables corresponding to 
different metallicities. 

2.2.3. Energy Generation 

Energy generation rates (e, in erg s" 1 g" 1 ) are computed by a 
custom nuclear reaction network. The network stores the rela- 
tive populations of the relevant nuclear species and derives from 
these and the reaction rates the energy generation rates from 
proton-proton, CNO, and He burning at supplied values of the 
primary and secondary variables. The 24 isotopic abundances 
appearing in Table 1 are followed individually in a fully im- 
plicit, time-dependent network modeled after the formalism in 
Hix & Thielemann (1999). Table 1 also lists the included reac- 
tions, for which the reaction rates are drawn from the canonical 
Caughlan & Fowler (1988) compilation. 

The nuclear reaction network is coupled to the stellar struc- 
ture equations through the secondary variable e and the abun- 
dances of the chemical species. These are stored in a separate 
array and the hydrogen, helium, and metal mass fractions, X, Y, 
and Z, respectively, are computed from this array, with scaling 
to the total density for use in the EOS and opacity routines. 

2.2.4. Convection and Mixing 



3 



Where the temperature gradient is too steep for radiative en- 
ergy transport, convective motions must carry the energy pro- 
duced in the stellar core. We assume that the isotopic abun- 
dances are fully mixed in convective regions and that convec- 
tive regions overshoot their formally defined boundaries. We 
follow the common practice (Maeder & Meynet 1987; Schaller 
et al. 1992) that assumes an overshooting distance of 0.20H p , 
where H p is the pressure scale height evaluated at the convec- 
tive boundary: 



H P =P 



dR 
dP 



(GM r /r 2 )p 



(5) 



Different assumptions for the overshooting distance may con- 
tribute modestly to the systematic uncertainty in the behavior 
of the stars in the HR diagram and their emergent spectra. 

2.2.5. Boundary Conditions 

Four boundary conditions are required to complete the sys- 
tem of four equations (Eq. 1 - 4). The outer boundary condi- 
tions constrain the temperature and pressure at the photosphere 
of the star. We assume that the temperature at the photosphere is 
the effective temperature given by the relation L = AirR^a^ T* n . 

To calculate the atmospheric pressure, we call a function that 
derives P atm = P(T,p,X,Y,Z) for the stellar photosphere. Near 
the stellar surface, we neglect changes in mass, radius, and lu- 
minosity and assume that M r = Af#, R = and L r = L*. This 
region obeys a simple relationship common to all stratified at- 
mospheres: dr = -npdr. Using this relation and combining 
equations 1 and 3, we derive a differential equation relating op- 
tical depth t and pressure P: 



dP GM* 



dr kRI 



(6) 



Then, using the Eddington approximation to radiative transport, 
we write the relation: 
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The boundary condition for the differential equation (8) is 
that the gas pressure vanishes at r = 0, leaving P = P a a = 
|(4(TsB/c)r 4 . At the photosphere, T = T e g and solving for t 
in equation (7) gives r = 2/3. We then integrate equation (6) 
from t = to t = 2/3, and find the value of P corresponding to 
T e ff at the photosphere. This calculation is done at each iteration 
in the approach to convergence. 

At the inner boundary we use the so-called "central expan- 
sions" (Kippenhahn et al. 1967), which fix the luminosity and 
radius at the inner mass zone. 

2.2.6. Initial Conditions 

The starting model for an evolutionary sequence is calculated 
at the zero-age main sequence (ZAMS). A mass and chemi- 
cal composition are specified, and a static structure model is 
calculated at t = 0. In this case, the starting models were the 
ZAMS models presented by TS. We use the initial composition 
X = 0.76, Y = 0.24 to represent the primordial composition. 
For illustration, we note that a decrease in X leads to system- 
atically higher core temperatures for stars powered by pp and 
CNO burning (e pp oc X 2 and ecNO oc XZq). 

For M ^ 15 Mq, ZAMS Pop III stars are not strictly metal- 
free except perhaps in the early stages of their formation (El 



Eid et al. 1983; Castellani et al. 1983; Marigo et al. 2000). 
The stellar cores are hot enough to accumulate a small abun- 
dance of 12 C via the triple-a process before and during main- 
sequence H-burning. Thus, the idea of a high-mass, metal-free, 
main-sequence star powered solely by pp burning is mislead- 
ing. In fact, these stars derive most of their energy from CNO- 
burning of H (see § 3.2 below). We must carefully follow the 
production of 12 C in the core to derive the initial conditions for 
a main-sequence (MS) Pop III model. In contrast, TS used ad 
hoc assumptions for the pre-MS synthesis of carbon, based on 
the earlier calculations of El Eid, Fricke, & Ober (1983). 

We start with the strictly metal-free ZAMS model from TS 
and let it produce 12 C in the core until CNO burning dominates 
the energy production. Starting from the Z = model, the star 
will evolve to slightly lower luminosity and r e ff until energy 
production from CNO burning exceeds that of the pp chains. 
At this point, the star begins to evolve to higher luminosity but 
still declines in T e g. The 12 C mass fraction at this transition 
is Zc ~ 10~ 12 (where Z ~ 0.02) and increases slightly with 
mass. We take the first point where dL/dt > to be the starting 
point of the evolutionary tracks. For M <; 15 M©, luminosity 
increases immediately, starting at Z = 0. 

3. RESULTS FOR STELLAR EVOLUTION AT ZERO METALLICITY 
3.1. Evolutionary Tracks 

We have calculated evolving stellar models for M = 1-100 
Mq from the ZAMS past the onset of helium burning. These 
tracks appear in Figure 1, which shows the tracks compared 
with a similar Pop II (Z= 0.001) zero-age main sequence forM 
= 1-100 Mq from Schaller et al. (1992). The large gain in 
effective temperature at Z = is readily apparent in this figure. 

We have done a simple comparison of our tracks to the sim- 
ilar work of Marigo et al. (2000) and Siess et al. (2002). Our 
tracks match theirs on the H-burning main-sequence, and for 
He-burning at M > 10 Mq. We find effective temperatures 
that are 3-5% lower and H-burning lifetimes that are 1-12 
% shorter (Table 2). These differences illustrate that individ- 
ual choices for parameters, numerical techniques, and the spe- 
cific implementation of physical processes (e.g., convection) or 
boundary conditions can create small differences in the final 
results. Our models diverge from Marigo et al. (2000) in the 
details of evolution in the HR diagram during the He-burning 
phases for M < 10 Mq, perhaps owing to differences in our 
treatment of convection. Since we are concerned with the spec- 
tral evolution for the first 10 Myr, which is dominated by the 
H-burning MS of massive stars, we have not attempted to re- 
solve these differences here. Simple testing indicates that the 
time evolution of ionizing photons and broadband spectra are 
similar between the two sets of tracks, except for a roughly 30% 
smaller He II ionizing photon production from our models, due 
to the slightly cooler effective temperatures and the sensitivity 
of 4 Ryd radiation to T e ff. We note also the studies of Pop III 
stars over a broad mass range by Schaerer (2002), and Bromm 
et al. (2001), who published static models of zero-age Pop III 
stars from 100- 1000 Mq. 

Our tracks do not complete core He burning because of a pe- 
culiar evolutionary feature noted also by Marigo et al. (2000). 
During core He burning, the core begins to expand and at some 
point reaches the boundary of the H-depleted region. At this 
point, the core H abundance rises from to ~ 10~ 4 from one 
discrete timestep to the next, owing to the assumption of instan- 
taneous complete mixing in the convective core. This assump- 
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tion is clearly inadequate here, since the added H is probably 
burned quickly at the edge of the convective region in a contin- 
uously evolving star. An ad hoc treatment of this phenomenon 
that forbids expansion of the convective core past this bound- 
ary is unlikely to represent accurately the true behavior of the 
stellar interior in this stage. Since this evolutionary phase rep- 
resents a small fraction of He burning, which is itself <;10% of 
the stellar lifetime, we terminate our tracks when the boundary 
of the convective core reaches the H-rich region. 

Figure 1 shows that, for M J> 12 Mq, the stars evolve uni- 
formly towards lower T e ^ and higher luminosity. These stars 
are powered by core CNO burning, catalyzed by the small C 
fraction built up there by 3a burning. For M <; 12 M Q , the stars 
start on the ZAMS with Z = and first evolve to higher T e ff via 
pp burning while the core C abundance rises. When the 12 C 
mass fraction reaches Z c ~ 10" 9 , the star is powered primarily 
by CNO burning and begins to evolve to lower T eS . 

To assess the differences in ionizing photon emission and 
some observational characteristics between Pop III and Pop II, 
we have constructed zero-age main sequences for initial car- 
bon abundances of Zq = 10~ 8 , 10" 7 , and 10" 6 . These appear in 
Figure 1 . These unevolved populations may represent the in- 
termediate stellar generations between Pop III and Pop II. We 
construct model spectra for these populations as for the Pop III 
evolving cluster (see § 4). 

We have neglected mass loss in calculating our stellar evolu- 
tion models and spectra. Mass loss is believed to be driven in 
massive stars by radiation pressure on blanketed metal lines in 
the outer reaches of the stellar atmospheres. In the metal-free 
case, this pressure is likely to be minimal (Kudritzki 2000), and 
the constant-mass assumption is a good approximation. Marigo 
et al. (2000) used a simple test of Eddington luminosity to argue 
that massive Pop III stars may experience modest mass loss late 
in their He-burning phases. Even if this occurs, it is unlikely to 
affect the composite spectra we produce here and can be safely 
ignored. The uncertain role of mass loss can be viewed as one 
of the major systematic uncertainties remaining in the study of 
metal-free stars. If it is strong, main-sequence mass loss drives 
the star to hotter T e ff on the main sequence - these effects can 
be seen in the tracks of El Eid et al. (1983) and Klapp (1983). 

3.2. Why Pop III Stars are Hot 

The distinctive spectral appearance of Pop III stars is due 
directly to their unusually hot surfaces (TS). These stars are 
roughly two times hotter and five times smaller than their Pop 
II counterparts of the same mass. There are two main reasons 
for these differences, both related directly to the chemical com- 
position. First, nuclear burning in a metal-poor isotopic mix 
must occur at higher temperatures to derive the same energy 
per unit mass as at higher metallicity. The change in nuclear 
energy production rates (e) is illustrated in Figure 2, with some 
characteristic parameters. The highest curve (green) is ecno for 
a metallicity of Zq = 10~ 6 . At this metallicity, CNO burning at 
a typical core temperature of 4 x 10 7 K can achieve the typical 
e for a particular location in the stellar interior. If this star is 
restricted to pp burning (black curve) or to CNO burning at a 
very low C abundance (red curve), it must do so at temperatures 
above 10 8 K. These latter two cases illustrate Pop III nuclear 
burning. This effect accounts for most of the large changes in 
stellar structure as metallicity approaches Z = 0. Second, and 
less important to the structure of Pop III stars, the opacity of 
stellar material is reduced at low metallicity, permitting steeper 
temperature gradients and more compact configurations at the 



same mass (see Figure 3). 

The effects of reduced opacity can be seen in Figure 3, which 
shows the run of six important physical quantities in Pop II 
(dashed) and Pop III stars (solid) of 15 Mq. Pop III stars are 
more centrally condensed, with hotter, denser cores, and hotter 
surfaces. Throughout 99% of the stellar material, the opacity 
(lower left) has a weak dependence on metallicity because it is 
dominated by scattering off electrons donated primarily by H. 
In the outer envelope, which comprises 1 % of the stellar mass, 
the opacity is provided by the bound-free and bound-bound 
transitions of H I, He I, and He II. Here the opacity differs sub- 
stantially between Pop II and Pop III, but it has little effect on 
the overall stellar structure. 

4. MODEL ATMOSPHERES AND SYNTHETIC SPECTRA 

4.1. Model Atmospheres 

The model atmosphere grid was calculated using the publicly 
available TLUSTY code (Hubeny & Lanz 1995). These mod- 
els are parameterized by the effective temperature at the stellar 
photosphere (r e ff) and the surface gravity (g). Model spectra at 
the grid points are calculated for A = 1 - 10, 000 A by the com- 
panion program S YNPLOT. For high accuracy and coverage we 
chose Ar e ff = 500 K and A(logg) = 0.05. Figure 4 shows the 
range of the model grid with the evolutionary tracks overlaid. 
At low temperatures (r e ff < 10,000 K) the model atmospheres 
are allowed to have small convective regions to promote con- 
vergence. This modification has little effect on their emergent 
spectra. Some stars develop convective envelopes during late H 
burning, but these regions lie far from the nuclear burning core 
and thus always maintain their metal-free composition. 

4.2. Population Synthesis 

To produce model spectra and calculate the cosmological ef- 
fects of Pop III stars we must create synthetic clusters from 
the discrete set of tracks presented above. We desire smoothly 
varying quantities such as ionizing photon flux and broadband 
colors, so we must convert the tracks at discrete masses into 
more continuous variables. 

Constructing a complete synthetic population from evolu- 
tionary tracks requires careful calculation of the tracks for 
masses that are not explicitly computed. For example, if one 
has tracks at 100, 80, and 60 Mq, extrapolating the behavior 
of the 90 Mq track from 60 and 80 Mq leads to large errors, 
and interpolation of the 100 and 80 Mq tracks is not accurate 
after the last timestep on the more massive track. We cannot 
simply use the explicitly calculated tracks to represent broad 
mass bins (AM ~20 Mq), because this leads to sharp discon- 
tinuities in the calculated spectra and derived quantities. We 
use small mass bins of 1 Mq, which are preferred for adequate 
accuracy and smooth behavior of the results. To resolve these 
issues we use a "lifetime-proportional" interpolation scheme to 
create synthetic evolutionary tracks for the intermediate masses 
between the explicitly calculated tracks. 

The evolutionary tracks at the explicitly calculated masses 
are followed beyond the onset of helium burning, a phe- 
nomenon that occurs at all masses and occupies < 10% of the 
stellar lifetime. Since later stages (C-burning, Si burning) are 
even shorter and contribute little to the spectral evolution of 
a cluster, we assume that He burning is the last evolutionary 
phase and that its end terminates the star's life. We also rely on 
the fact that the changes in shape and time-dependence of the 
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tracks vary smoothly with mass. We can use these important 
features to construct interpolated tracks at intermediate masses. 

Taking the run of T e s(t) and L(t ) to be a track, we first remap 
the explicit tracks to a new time coordinate - the fraction of 
the stellar lifetime achieved to that point, ranging from to 1 . 
Then, a vector of lifetimes for each mass bin is interpolated 
from the lifetimes of explicit tracks. Finally, we interpolate 
intermediate tracks between the explicit tracks at the same frac- 
tion of their lifetimes. For instance, to find (L, T e s) for a 90 M Q 
star at 70% of its interpolated lifetime, we interpolate in mass 
between the 80 M and 100 M Q explicit tracks at 70% of their 
lifetimes. This scheme minimizes interpolation errors and gen- 
erates smoothly varying model spectra and derived quantities. 

Using the interpolated tracks and model atmospheres, we 
construct synthetic stellar populations by specifying a stel- 
lar initial mass function (IMF) and adding stellar spectra to- 
gether. Because we lack detailed knowledge of the IMF for 
Z = star formation, we assume that the IMF is a power law, 
dN/dM oc M~ a , with the Salpeter slope, a = 2.35. To populate 
the IMF, we construct 99 mass bins and populate them with 
the number of stars of that mass, normalized to a total mass 
of 10 6 Mq, which is standard practice. The mass bins are 1 
Mq wide and range from 1-100 Mq. Model spectra from the 
atmosphere grid are assigned to the bins and summed over all 
masses. The product of this process is a composite spectrum 
for the whole cluster. Schaerer (2002) and Bromm et al. (2001) 
performed similar calculations for top-heavy IMFs ranging up 
to 1000 M Q . 

We adopt a model of Pop II starburst clusters as a reference 
spectrum for comparison with our evolving Pop III models. In 
the present paper, this reference model represents a baseline 
spectrum for comparison to the Z = models. In Paper II, we 
use this spectrum as a proxy for the "second generation" of stars 
incorporating the metals produced by the first. For the refer- 
ence population, we use a model with Z = 0.001, M = 10 6 M 
in stars, and a Salpeter IMF, that was custom-generated by the 
Starburst99 population synthesis code (version 3.1, Leitherer 
et al. 1999). Because this suite of models is widely used in 
the community, we compare all our model spectra and obser- 
vational predictions to a Pop II model generated by their code. 
We note further that we do not include Wolf-Rayet (WR) stars 
in these comparison models to make our comparisons as direct 
as possible. WR stars are known to produce large amounts of 
He II ionizing photons in the local universe (Schaerer & Vacca 
1998), but their at high redshift (z > 5) and low metallicity is 
regarded as unlikely since the WR phenomenon requires metal- 
line opacity to drive strong stellar mass loss. Recent studies of 
stellar mass loss at low metallicity (Kudritzki 2002) have shown 
that its importance is minimized at extremely low metallicity, 
with a factor of 10 drop in mass-loss rate from Z = 10" 3 to 10" 4 . 
Nonetheless, we include in our comparisons below some results 
for Pop II with WR stars. 

For modeling cosmological reionization, radiative feedback, 
and the observational signatures of metal-free stars, we make 
use of the important quantities g,, which express the ionizing 
photon output per second in the H I, He I, and He II continua 
(i = 0, 1, 2, and hv t = 13.6, 24.6, 54.4 eV, respectively): 

f°° F 

Qi=4nRl / ^-dv, (8) 
Jv, hv 

where F v is the spectrum of the star in units of erg cm" 2 s" 1 
Hz" 1 . The time-dependent values of go an d Qi are displayed 



for the individual stars in Figure 5 to show the behavior with 
mass. 

Once the individual g, have been weighted by the IMF, we 
arrive at the cluster g, seen in Figure 6. There are several fea- 
tures of this figure that deserve comment. First, note the starting 
(ZAMS) values of go and g 2 for the Pop III and comparison 
Pop II clusters. The Pop III cluster has go that is 60% larger 
than the Pop II cluster. Second, the ionization of the Pop III 
cluster persists at a high level far longer than the Pop II cluster, 
owing to the strong gain in ionization forM = 10-25 Mq stars. 
The Pop III cluster takes 10 Myr to decline to 1/10 its ZAMS 
value in go, twice as long as Pop II. Finally, note the strong ion- 
ization above the He II edge at 4 Ryd. At the Pop III ZAMS, we 
find that QiIQq = 0.02, with a lifetime of 2.5 Myr, characteristic 
ofM>40-50M Q stars (see Figure 5). This form of the output 
from our models can be placed directly into cosmological cal- 
culations of IGM reionization that incorporate realistic models 
of evolving spectra from Pop III stars (see Paper II). 

For the ZAMS spectra at Z c = 10" 8 , 10" 7 , and 10" 6 , we 
find go = (8.9,8.4,7.9) x 10 52 s" 1 , respectively, intermediate 
between the Pop II and Pop III values. There is a sharp de- 
crease in their g 2 , however, with g 2 = 8.4 x 10 49 , 2.1 x 10 48 , 
and 4.8 x 10 47 s" 1 , respectively. These quantities are compared 
with the evolving g 2 for Pop III and Pop II in Figure 6. 

In Paper II, we examine the effects of Pop III clusters on the 
reionization of the IGM. These tests use a quantity r\, which ex- 
presses the ratio of He II to H I column densities in an ionized 
IGM (Fardal, Giroux, & Shull 1998). This ratio can be pre- 
dicted for a specific distribution of IGM clouds, given an input 
source spectrum, which is assigned an intrinsic rj (prior to IGM 
filtering) based on the shape of its spectrum above 1 Ryd. A 
composite QSO (Zheng et al. 1997; Telfer et al. 2002) with 
a power law spectrum, F v ~ v~ a , and a = 1.8 has an intrinsic 
r\ = 20. At the ZAMS, the Pop III cluster has 77 = 10, harder 
than the mean QSO and similar to a power-law spectrum with 
a = 1.3. For the low-metallicity ZAMS at Z = 10" 8 , 10" 7 , and 
10" 6 , we find intrinsic r\ = 33, 80, and 180, respectively. We 
make use of this ratio in Paper II to explore the fate of IGM 
regions ionized by Pop III stars. 

5. OBSERVATIONAL SIGNATURES 

5.1. General Approach 

Observational searches for high-redshift galaxies work best 
when equipped with diagnostic tools that provide detailed rela- 
tionships between observable signatures and underlying phys- 
ical phenomena. To this end, we use our models of Pop III 
clusters to derive predictions for the observed spectral energy 
distributions (SED) of the first stars. We consider two classes 
of observation; the first seeks the nebular emission line signa- 
tures of Pop III stars, and the second exploits their broadband 
SED and colors. We take up the nebular emission first because 
its presence affects the broadband signatures as well. 

5.2. Emission-line Diagnostics 

Emission-line signatures are proving to be the primary means 
of finding galaxies at redshifts beyond z = 5. Deep, blind 
searches for Lya emitters have turned up large populations of 
galaxies at these epochs. These surveys use either blind multi- 
slit spectrographic exposures to detect single emission lines 
over a broad redshift range (Hu et al. 1999), or they take nar- 
rowband images to probe deeply at selected redshifts (Rhoads 
et al. 2000). 
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TGS derived a relationship for the observed He II emission- 
line flux from metal-free stars as a function of star formation 
rate (SFR). With the evolved spectra of metal-free stars we can 
update these earlier calculations. We assume that few He II ion- 
izing photons escape the galaxies where they are produced (/ esc 
= 0.025; see Paper II) and that the sources themselves contain 
no dust. We assume Case B recombination at 20,000 K, corre- 
sponding to a higher nebular temperature in low-metallicity gas. 
As in TGS, we adopt j^ml 'ha = 0.66, jmo/ha = 4.6, and 
73203/ 7Ha = 0.28 (Seaton 1978). TGS used an ad hoc scale fac- 
tor, /evoi, to account for stellar evolution effects, which reduce 
the He II ionizing flux over time as high-mass stars evolve off 
the main sequence. We need not make this approximation here, 
since we have directly calculated evolving spectra of metal-free 
stars. For both continuous and instantaneous star formation, 
we assume that L Ha = 1-27 x 10~ 12 go erg s" 1 (i.e., there are 
0.424 Ha photons emitted for each H-ionizing photon). Simi- 
lar coefficients are derived for the other H I lines from Case B 
recombination theory (Osterbrock 1989). Note that Leitherer et 
al. (1999) use a coefficient of 1.36 x 10~ 12 , corresponding to 
nebular gas at 10 4 K. We then multiply the Ha luminosity by 
the emissivity ratio _/,/ j Ha to derive L, in a He II line. In Fig- 
ure 7 we plot the predicted line fluxes for Lya and He II A 1640 
as a function of redshift. 

For clusters with significant ionizing photon production, we 
must also account for continuous nebular emission from the H I 
two-photon process and direct recombinations to upper levels 
of H I, He I, and He II. We construct a wavelength-dependent 
emission spectrum, in units erg s" 1 A" 1 : 



rcont 



>i n (1 - ./esc)2o, 



(9) 



where the term 7,, in units erg cm 3 s 1 Hz 1 , incorporates con- 
tributions from H I, He I, He II, and the H I two-photon contin- 
uum and aL^ is the Case B recombination rate coefficient for 
hydrogen. The emission coefficients for an electron tempera- 
ture of 20,000 K are taken from Aller (1984, Table 4.9). This 
additional emission has a large effect on the broadband spec- 
tral energy distribution of the synthetic clusters, shifting en- 
ergy from the Lyman continuum to the rest-frame optical bands. 
This effect is included in the broadband color and equivalent 
width predictions presented below. 

In general, we can parameterize the escape of radiation with 
the quantity / esc , the fraction of ionizing photons that escape 
from the host halo into the IGM. To be consistent with the typi- 
cal reionization models presented in Paper II, we uniformly as- 
sume / esc = 0.05 for 1 - 4 Ryd and / esc = 0.025 for hv > 4 Ryd 
(Dove, Shull, & Ferrara 2000; Tumlinson et al. 1999), leaving 
most ionizing photons at all wavelengths to excite observable 
nebular emission. We return to the issue of the escape fraction 
in Paper II, where we discuss justifications for different escape 
fractions at 1 and 4 Ryd. 

In Figure 7 we show the flux of the A 1640 and Lya lines from 
the Pop III cluster. The green curves represent continuous star 
formation at 1 and 40 M Q yr" 1 . The higher rate corresponds to 
the star formation rate (SFR) inferred by Hu et al. (2002) for a 
Lya emitter discovered at z = 6.56 (Lya flux marked with the 
green square). The red curves represent the He II A1640 A flux 
from an instantaneous burst of M = 10 6 M at t = (solid) and 

4 We note here that our Lya predictions lie well below those of Schaerer (2002) for similar conditions. We use an average flux in the 1200 - 1250 A range to 
represent the continuum flux. This method approximates the observed quantities in surveys such as LALA. If intrinsic Lya absorption is present, our prediction is an 
overestimate. We can reproduce Schaerer's prediction only by assuming the continuum flux is represented by the core of the model Lya absorption line. 



1 Myr after the burst (dashed). Figure 7 also shows the imaging 
and spectroscopic sensitivity requirements for the planned 6 m 
aperture Next Generation Space Telescope (NGST). These as- 
sumed limits are from the NGST design specifications and are 

2 nJy for low-resolution (R = 5) imaging and 100 nJy for spec- 
troscopy (R = 1000), and both assume S/N = 10 and a 100,000 
s exposure time. These comparisons show that the planned 
NGST will easily discover this signature of the first stars. 

In Case B recombination theory, roughly two thirds of all H 
recombinations result in the emission of a Lya photon (Spitzer 
1978), and these are known to escape from galaxies and serve 
as a strong signature of high-redshift star formation (Hu et al. 
1999; Rhoads et al. 2000). For the Pop III cluster presented 
here, the Lya line is intrinsically stronger than in Pop II. Fig- 
ure 8 shows the predictions for the time-dependent Lya equiv- 
alent width, WLy Q , for Pop II and Pop III assuming Case B re- 
combination theory and T = 20000 K. These curves are tech- 
nically upper limits since the Pop II cluster may suffer dust 
attenuation and the Pop III Lya line may be absorbed by neu- 
tral H in the pre-reionization IGM (cf. Loeb & Rybicki 1999). 
Equivalent width has the additional benefit of being indepen- 
dent of the total cluster mass for an instantaneous burst. Thus, 
we confirm the suggestion of other studies (TS, Bromm et al. 
2001; Schaerer 2002) that the Lya signatures of the first stars 
will be unusually strong 4 . The Large Area Lya Survey (LALA; 
Malhotra & Rhoads 2002) found a median W Lya of 400 A in a 
sample of 157 galaxies with z = 4.5. More than 60% of their 
detections have Wi, y a > 240 A, the approximate maximum at- 
tainable by Pop II clusters (Chariot & Fall 1993). 

We also plot in Figure 8 the instantaneous equivalent widths 
W Lya = 520, 430, and 360 A, for the Z c = 10" 8 , 10- 7 ,andl0" 6 
ZAMS spectra, respectively. The points forZ c = 10" 7 and 10" 6 
straddle the median LALA equivalent width, indicating that this 
sample of Lya emitters may have a low metallicity. These re- 
sults suggest that populations of metal-free stars, or perhaps 
stars of the first metal-enriched generation, have already been 
found. In future ground- and space-based searches beyond the 
optical bands, Lya emission will probably serve as the "sign- 
post" of metal-free galaxies. This result could then be con- 
firmed through follow-up detection of the He II lines, which 
are weaker than Lya and unique to the metal-free regime. 

The sensitive dependence of the Lya equivalent width on 
metallicity (and IMF; see Bromm et al. 2001) suggests that 
the metallicity evolution of star formation may appear in the 
observed Lya lines at high redshift. If there is a well-defined 
Pop III epoch, characterized by a rapid transition to Zq = 10~ 6 
or higher as the second stellar generation forms from gas en- 
riched by the first, there could be a bimodal distribution of Lya 
equivalent widths with redshift. The peak may appear near the 
~ 1000 A prediction here, nearer the <~ 3000 A predicted for 
the high-mass IMF adopted by Bromm et al. (2001), or some- 
where in between. Thus, Lya itself may serve as a diagnostic 
of the evolution of metals in galaxies, up to at least the redshift 
of H reionization, where the damping wing of Gunn-Peterson 
absorption in the IGM may play an uncertain role in obscuring 
(Miralda-Escude 1998) or transmitting (Loeb & Rybicki 1999; 
Haiman 2002) the galactic Lya line. 

5.3. Broadband Spectra and Colors 
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Broadband filter search techniques have proven successful in 
finding high-redshift galaxies. In particular, the "Lyman-break" 
technique exploits the fact that the IGM absorbs the spectrum 
of a galaxy shortward of 1216(l+z) A. This powerful technique 
has uncovered thousands of galaxies at z > 3 (Steidel et al. 
1999). 

The high effective temperatures of Pop III stars make their 
integrated cluster spectra intrinsically bluer than their Pop II 
counterparts. However, the most distinctive feature of the Pop 
III spectrum is the large gain in H I and He II ionizing photons, 
which unfortunately are emitted in the spectral region that is 
heavily absorbed by the IGM. To account for this obscuration, 
we adopt a stochastic model of the IGM based on the work of 
Fardal, Giroux, & Shull (1998), who used Monte Carlo tech- 
niques to construct model sightlines containing realistic dis- 
tributions of H I Lya forest, Lyman limit, and damped Lya 
clouds. From a sample of 1000 sightlines we construct a mean 
transmission function for three redshifts, z = 4, 5, 6. This real- 
istic treatment of the IGM leaves approximately ~ 30% of the 
emergent flux between 912(l+z) and 1216(l+z) A. This mean 
transmission function is included in the broadband signatures 
derived here. 

In Figure 9, we illustrate Pop II and Pop III spectra at high 
redshift. We also plot the ranges of the standard photometric 
bandpasses. Emission lines are included in these colors and 
make a modest contribution. Although the Pop III spectrum 
is intrinsically bluer in the accessible spectral regions, the re- 
processing of ionizing photons into nebular emission makes up 
for this difference. Figure 10 illustrates this similarity between 
Pop III and Pop II over a range of redshifts and cluster ages (f 
= 0-10 Myr). This substantial degeneracy in the broadband 
colors may make distinguishing Pop III more challenging for 
current and planned optical and IR instruments than was pre- 
viously appreciated. Confusing color signatures may leave the 
unusual Pop III emission-line signatures as the only means of 
identifying the first stars. 

6. CONCLUSIONS 

We have presented evolving spectra of metal-free stellar pop- 
ulations, based on newly calculated evolutionary tracks. We 
find that the evolution of Pop III stars follows the general pat- 
terns obeyed at higher metallicity, but with an overall shift to 
higher temperature. This gain in core and surface temperatures 
is primarily a result of the restricted abundance of 12 C in pri- 
mordial stars. We have used a grid of non-LTE model atmo- 
spheres to produce evolving spectra of synthetic Pop III clus- 
ters. These models have been used to compute the broadband 
colors and emission-line signatures from metal-free stellar clus- 
ters at high redshift. 

Our specific conclusions fall into two categories. First, the 
unique composition of the first stars has the following impor- 
tant effects on their spectra: 

• They produce 60% more H I ionizing photons than their 
Pop II counterparts. 

• They can produce up to 10 5 times more He II ionizing 



photons than Pop II, which may lead to unusual radiative 
feedback effects on the IGM. 

We wish to reemphasize here the uncertainty associated with 
stellar mass loss in evaluating the latter conclusion. If the "sec- 
ond generation" of stars has Z ~ 0.0001-0.001 and perhaps 
significant populations of WR stars, the gain in He II ionizing 
photon production for a continuous star formation model in- 
cluding WR stars is reduced to about a factor of 10 (see Figure 
6). However, since the importance of WR stars at extremely 
low metallicities is in doubt, we conclude that Pop III stars will 
be far more efficient than their Pop II counterparts at ionizing 
He II. The consequences of our models for reionization are ex- 
plored more fully in Paper II, which uses these model spectra 
and a semi-analytic reionization model to assess the importance 
of the first stars for full or partial reionization of H I and He II 
in the IGM. 

The peculiar properties of metal-free stars entail unusual ob- 
servational signatures: 

• They are expected to have unusually high equivalent 
widths of Lya and the He II recombination lines, the 
most distinctive signature of Pop III stars and their neb- 
ular emission. These lines will probably be the primary 
means of detecting and identifying the first stars. 

• While their intrinsic spectra are significantly bluer than 
their Pop II counterparts, their broadband colors are sim- 
ilar, owing to the reprocessing of Lyman-continuum pho- 
tons into continuous nebular emission. 

Recent Lya surveys suggest that metal-free stars may have 
already been found. In future ground- and space-based searches 
beyond the optical bands, Lya emission will probably serve as 
the "signpost" of metal-free galaxies. Follow-up detections of 
the He II lines could confirm this result. The substantial degen- 
eracy in the broadband colors will probably leave emission line 
techniques as the best way of distinguishing the first stars. 

We are grateful to Mark Giroux for calculating the 
IGM transfer functions presented in Section 4, and to 
Claus Leitherer and his collaborators for the public dis- 
tribution of their Starburst99 population synthesis code 
(http://www.stsci.edu/starburst99). Comments from John 
Stocke, Nick Gnedin, and Andrew Hamilton improved the 
manuscript. We happily acknowledge the anonymous referee 
who prompted us to clarify the discrepancy between our emis- 
sion line predictions and those of Schaerer (2002), and who 
provided other useful comments. The NGST design speci- 
fications were taken from the NGST website in April 2002 
(http://www.stsci.edu/science/sensitivity). The OPAL project 
at Lawrence Livermore National Laboratory is gratefully ac- 
knowledged for their opacity data and interpolation codes 
(http://www-phys.llnl.gov/Research/OPAL/opal.html). This 
work was partially supported by an astrophysical theory grant 
to the University of Colorado by NASA (NAG5-7262). 
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Table 1 

Reactions included in the stellar evolution models. 

Class Reaction 
~PPl l H+ l H^ z H+v e + e + 

2 H+ 1 H^ 3 He + 7 
3 He + 3 He 4 He + 2 ! H 

PPII 3 He + 4 He -> 7 Be + 7 

7 Be + e" — > 7 Li + v e 
7 Li+ 'H^2 4 He 

PPIIII 7 Be + 'H 8 B + 7 
8 B -> 8 Be* + e + + v e 
8 Be* -> 2 4 He 

CNO 12 C + 'H -> 13 N + 7 

13 N -> 13 C + e + + v e 
13 C+ 'H-^ 14 N + 7 
14 N+ >H-> 15 + 7 
15 -> 15 N + e + + v e 
15 N+ >H-> 12 C + 4 He 
15 N+ >H-> ls O + 7 
16 0+ 'H^ 17 F + 7 
17 F -> n O + e + + v e 
n O+ 'H^ 14 N + 4 He 
17 0+ 'H^ 18 + e + + v e 

3a 4 He + 4 He -> 8 Be* 

8 Be* -» 2 4 He 
8 Be* + 4 He 12 C + 7 

a capture 12 C + 4 He -> lfi O + 7 
lfi O + 4 He 20 Ne + 7 
20 Ne + 4 He -» 24 Mg + 7 
14 N + 4 He -> 18 F + 7 
18 + 4 He -> 22 Ne + 7 
22 Ne + 4 He -> 25 Mg + n 



Table 2 
pop iii h-burning lifetimes 



Mass (M Q ) 


th (yr) 


1 


5.87xlO y 


3 


2.01 xlO 8 


5 


6.03 xlO 7 


7 


3.11xl0 7 


10 


1.76xl0 7 


12 


1.40xl0 7 


15 


l.llxlO 7 


20 


8.47xl0 6 


30 


5.86xl0 6 


50 


4.09 xlO 6 


70 


3.42xl0 6 


100 


2.93 xlO 6 
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5.0 4.8 4.6 4.4 4.2 4.0 3.8 

log T eff [K] 

FIG. 1 . — The Z = stellar evolution tracks from 3 - 100 Mg . The 1 Mq track has been omitted for clarity. These tracks terminate when the convective core 
expands to reach the H-rich region during core He burning. To illustrate the stark differences between Pop III and their metal-enriched counterparts, we mark in red 
the ZAMS for Pop II (Z = 0.001) from Schaller et al. (1992). The Pop III ZAMS is marked in blue, and the ZAMS for Z c = 10~ 8 , 10~ 7 , and ICr 6 are plotted with 
dashed lines, with metallicity increasing to the right. 



11 




0.00 0.05 0.10 0.15 0.20 

T/10 9 [K] 

FIG. 2. — Energy generation rates (e) illustrate why metal-free stars must attain unusually high temperatures to support their mass against gravity. The black curve 
shows the pp rate for a typical H mass fraction X = 0.7. The green and red curves show the CNO rate for a fiducial low metallicity and for Pop III, respectively. The 
horizontal dashed line marks a typical (but arbitrary) rate for the interior of a massive (M > 20 Mq ) star. The circles mark the points where the respective burning 
processes produce the desired total e. See text for discussion. 
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FIG. 3. — A comparison of the interior structure of a Pop III star (solid curves) with a Pop II star (dashed curves), both of M t = 15 Mq. The Pop III star is smaller, 
hotter, and denser. This comparison demonstrates that the gain in r e ff arises not from changes in the opacity with metallicity, which differs substantially only in the 
outer 1% of the total mass, but rather from changes in the available fuel source in the stellar core. See § 3.2 for discussion. 
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log t [yr] 

FIG. 5. — Ionizing photon production rates go (1 Ryd; upper panel) and Qi (4 Ryd; lower panel) for individual stars from the tracks in Figure 1. The tracks 
are labeled with their stellar mass in Mq. High-mass, zero-metallicity stars are strong sources of H I and He II ionization. In particular, they emit 10 5 times more 
radiation in the He II continuum than Pop II stars. 
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FIG. 6. — Upper left: composite spectra for Pop II (black) and Pop III (red) clusters (M = 10 6 Mq in stars) at 1 and 15 Myr. At 15 Myr the Pop II spectrum has 
faded in H I ionization, but the Pop III cluster is still a significant source of H I ionization owing to the presence of stars with M = 10—15 Mq and T e n > 50,000 
K. No nebular emission is included here. Lower left: example spectra for the continuous star formation case at 15 Myr. Upper and lower right: cluster Qo and 
Q2 corresponding to the instantaneous/continuous cases at left. H I ionization from the Pop III cluster is 60% stronger than Pop II. Pop III emits 10 5 times more 
He II ionizing photons than Pop II with the same IMF, total mass, and excluding WR stars. If WR stars are included, this gain in He II ionization is smaller but still 
substantial. This key result has potentially large effects on the IGM (see Paper II). In the upper right panel we mark with filled squares the instantaneous values of 
Q2 from the zero-age main sequences with Zq = 10~ 8 , 10~ 7 , and 10~ 6 (plotted at an arbitrary time). These points show the sharp decline in He II ionizing photon 
production when small abundances of 12 C are included. The right panels also show the Starburst99 results with WR stars (dotted lines), for comparison. See § 4.2 
for discussion. 
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Z 

FIG. 7. — Here we show the predicted fluxes for hya and He II A1640 from the Pop III cluster, hya emission may serve as "signposts" pointing the way to 
metal-free star formation, which could then be confirmed by the detection of He II lines. The green square marks the observed hya flux detected from the z = 6.56 
galaxy found by Hu et al. (2002), and the vertical arrow connects this point to the predicted He II flux curve if this galaxy harbors an continuous Pop III starburst 
with SFR = 40 Mq yr 1 . The lower green curve corresponds to SFR = 1 Mq yr 1 . The red curves mark the He II A 1640 flux for an instantaneous burst at and 1 
Myr (solid and dashed, respectively). The blue curves plot the imaging and spectroscopic sensitivities for NGST. At the top are wavelength ranges of the common 
photometric bandpasses. The horizontal lines extending to z = 5 represent the sensitivity limits of recent emission-line surveys of galaxies: from top to bottom, 
Rhoads et al. (1999), Hu et al. (1999; imaging data). Stern et al. (1999), and Hu et al. (1999; spectroscopic data). 
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FIG. 8. — Predicted equivalent width, W\, in the hya emission line for the Pop III synthetic cluster (a Salpeter IMF up to 100 Mq). The Pop III equivalent width 
is unusually large, owing to the harder Pop III spectra. This horizontal dashed line represents the 80 A detection limit of the Large Area hya survey (Malhotra & 
Rhoads 2002). The horizontal dotted line marks 400 A, the median equivalent width for their sample of 157 galaxies. The filled squares mark the instantaneous 
Wiyty from the low-metallicity ZAMS, placed arbitrarily at t = 10 yr for clarity. Pop III clusters may be identified by their large Lya equivalent widths before their 
continuum spectra or broadband colors are measured. However, we note that at Pop II metallicities, a top-heavy IMF or one with stars more massive than 100 Mq 
can mimic this effect. 
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FIG. 9. — Example observed spectra, at time t = 0, for the Pop III (left) and Pop II (right) cluster at z = 5. Shortward of the rest-frame Lyman limit, the spectra 
are marked with long-dashed lines. The black lines mark the intrinsic stellar spectra, and the red lines include the nebular continuum emission. Labeled along the 
top are the common optical and infrared photometric bandpasses. The dotted line in the right panel repeats the Pop III spectrum at left for comparison. The spectral 
regions where Pop III luminosity exceeds that of Pop II lie in the range heavily attenuated by the IGM (A < 1216(1 +z) A), which is represented here by a mean 
transmission for 1000 sightlines simulated by a Monte Carlo model of the IGM (Fardal, Giroux, & Shull 1998). Note the absence of the He II recombination lines 
in the Pop II spectrum and the similarity between the observed spectra when nebular continuum emission is included. 
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FIG. 10. — The important broadband color diagnostics for the Pop III cluster (dotted lines and filled squares), compared to Pop II (solid lines and filled circles) 
from the ZAMS to 10 Myr in age, at redshifts z = 4, 5, and 6 (blue, green, and red, respectively). The tracks start at the filled symbols. Pop III clusters are intrinsically 
fainter and bluer in the infrared bands, but the influence of nebular continuum emission for t < 10 Myr hide this difference and give the Pop III and Pop II clusters 
similar color signatures. Thus, the recombination lines of H I (Lyo) and He II (A1640) may prove to be the best means of distinguishing Pop III stars. 



